32 Bin Near-Infrared Time-Multiplexing Detector with Attojoule Single-Shot Energy 
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We present two implementations of photon counting time-multiplexing detectors for near-infrared 
wavelengths, based on Peltier cooled InGaAs/InP avalanche photo diodes (APDs). A first imple- 
mentation is motivated by practical considerations using only commercially available components. 
It features 16 bins, pulse repetition rates of up to 22 kHz and a large range of applicable pulse widths 
of up to 100 ns. A second implementation is based on rapid gating detectors, permitting deadtimes 
below 10 ns. This allows one to realize a high dynamic-range 32 bin detector, able to process pulse 
repetition rates of up to 6 MHz for pulse width of up to 200 ps. Analysis of the detector response 
at 16.5% detection efficiency, reveals a single-shot energy resolution on the attojoule level. 



I. INTRODUCTION 

Conventional semiconductor-based single photon de- 
tectors, i.e. single-pixel Geiger-mode avalanche photo 
diodes (APDs) , are used in very different domains where 
high sensitivity is demanded. Naturally they are binary 
detectors, giving a click or no-click in response to an 
incoming light pulse. While this behavior is desired in 
fields like quantum information or quantum cryptogra- 
phy, it represents a restriction when light pulses contain- 
ing a few to hundreds or thousands of photons have to 
be measured. The restricted dynamic requires the use of 
an attenuator and numerous measurement repetitions in 
order to infer the average pulse energy p[|. Therefore a 
detector with the same sensitivity as single-pixel photon 
counting detectors, but with a larger response dynamic 
is desirable. It would strongly reduce the measurement 
time and allows one to infer single-shot estimates of the 
pulse energy. This can, in particular, be an advantage in 
biological imaging, where the number of applicable rep- 
etitions is restricted by the sensitivity of the specimen, 
or, in the characterization of mesoscopic quantum states 

i- 

Detectors that are capable of estimating the number of 
photons in a pulse on a single-shot basis are commonly 
called " photon number resolving detectors" . The resolu- 
tion on that number strongly depends on the efficiency 
of the detector. Examples for the visible wavelength 
regime are the silicon photo multiplier (SiPM) 0, Ij], 
time-multiplexing detector 0, 0j and visible light pho- 
ton counter (VLPC) 0, 1^ , and for the near-infrared the 
multi-pixel array frequency upconversion combined 
with SiPM [l3], avalanche height discrimination detec- 
tor [ll| and transition edge sensors (TES) jl^, [ll]- In 
the near-infrared, detectors either suffer from high dark 
count rates [oj, low efiiciencv flC | , p oor dynamic flT'l or 
depend on cryogenic cooling [12l [l3| . 

In this paper we present two implementations of time- 



multiplexing detectors for the near-infrared wavelength 
regime (1100 - 1650 nm), exhibiting low dark count, com- 
pact cooling, high dynamic and the capability of process- 
ing high pulse repetition rates. 

Time-multiplexing means that a pulse with a large 
number of photons is split into well separated weaker 
pulses by means of fiber delay loops (FiglT|), subsequently 
detected by single-pixel InGaAs/InP APDs. This kind of 
detectors has been implemented so far only in the visi- 
ble wavelength regime using Silicon APDs Up to 
three fiber loops (22, 44 and 88 m), distributing the initial 
photons on 16 different bins, were demonstrated. Due to 
high fiber loss at visible wavelengths (> 6 dB/km), these 
detectors are not efficiently scalable to a higher number 
of loops. 

Our first detector version uses commercially available 
Geiger-mode APDs (ID200, IDQ) and represents an easy- 
to-implement and easy-to-control scheme. The second 
version applies rapid gating detectors [l3 - [l7l | with dras- 
tically decreased afterpulsing, allowing system scalability 
due to low fiber loss at near-infrared wavelengths (< 0.2 
dB/km). 

The remainder of this paper is organized as follows: 
Section |ll] describes both setups and their characteriza- 
tion using coherent pulses. In Section iHll we evaluate the 
energy resolution for single and multi-shot measurements 
and section Hv] contains our conclusion. 
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II. SETUP 

The common setup for both multiplexing detectors is 
shown in Fig[T] It consists of two InGaAs/InP APDs, a 
time-multiplexer constructed by single-mode fibers and a 
control and acquisition system. The latter is responsible 
for the synchronized timing of pulses and gates, and the 
acquisition of detections. The multiplexing part includes 
m fiber loops (m=3 for the first and m=4 for the sec- 
ond implementation) fusion spliced via 50/50 couplers. 
The initial pulse is split into 2™+^ weaker pulses, subse- 
quently called bins, 2"^ directed to each detector. In the 
ideal case loop i has at least twice the length of the previ- 
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FIG. 1: Schematic of the common setup, comprising a time- 
multiplexer (m fiber loops), two InGaAs/InP detectors and a 
central control and acquisition unit, which is responsible for 
accurate synchronization of laser pulses and detector gates as 
well as for the registration of detections. 



loop lengths, however, have to be chosen in accordance 
with the detector characteristics. Here, the tolerable af- 
terpulsing is crucial. To keep its impact low, the distance 
between adjacent bins (deadtime) has to be chosen suffi- 
ciently large. 



A. Conventional gating detectors 

Two Peltier cooled InGaAs/InP APDs (ID200,IDQ), 
using a conventional gating technique with applicable 
gate widths between 1 ns and 100 ns, are gated exclu- 
sively during the arrival of each bin, see FigH) This is 
facilitated by an universal timing device (UTD) (Chrono- 
logic Legato), which sends a trigger signal to the laser 
and with a certain delay a gate pattern to the APDs. 
The timing of the pattern can be defined by software in 
correspondence to the bins of the multiplexer. The tim- 
ing accuracy of the module is typically 300 ps. The UTD 
also acquires the number of detections by means of two 
additional ports, subsequently stored on the computer 
(USB connection). 

We choose a deadtime of 5 /is at a detection efficiency 
of rj =10%. The afterpulse probability in a 20 ns gate 
after a detection gate is equal to 9%. The dark count 
probability is 8 • 10~^ per ns of gate. The deadtime set- 
ting demands a first loop of at least 1 km. In total we 
implement 3 loops of 1 km, 2 km and 5 km, yielding tem- 
poral delays of 5, 10 and 25 /is and 16 bins. The overall 
transmission loss, averaged over all bins, is 1.44 dB. 

We characterize the detector by sending coherent 
pulses of 1 ns-width provided by a diode laser at 1559 nm. 
They can be attenuated down to a well defined average 
number of photons per pulse by means of a calibrated 
variable attenuator. The maximally processable laser 
pulse repetition rate is 22 kHz, limited by the total length 
of the pulse train of 45 /is. 
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FIG. 2: Timing of conventional gating time-multiplexing de- 
tector. The universal timing device (UTD) sends laser trigger 
pulses with a fixed frequency, which generates 16 bins after 
the multiplexer. Gate signals are sent with an adequate delay 
in order to activate the detector at the pulse arrivals (bins). 
The number of occurred detections is then registered by the 
UTD. 



We acquire data typically for a few seconds until suf- 
ficient statistics are obtained. The laser repetition rate 
in these measurements is 10 kHz. Fig |3] shows the prob- 
ability distributions p(n\fi) (red circles with error bars) 
of the number of detections n obtained by sending weak 
coherent pulses of average photon numbers n — 7, 13, 50 
and 100, using a gate width of 20 ns. In order to evaluate 
our results, we perform a simple Monte Carlo simulation, 
modelling the fiber couplers as perfect 50/50 beam split- 
ters and accounting for the overall transmission loss in 
the multiplexer as well as dark counts and detector ef- 
ficiency. The /X— value that best fits the experimental 
results (plotted in FigI3]) is on average 8% higher than in 
the pulse that is sent on the detector. This offset can be 
explained by the fact that afterpulsing is not accounted 
for in the simulation. We conclude that even without a 
specific model for afterpulsing the detector behavior can 
be simulated in good agreement with the actual detector 
response. 

A possible extension of this detector, realizing a forth 
loop, would require an additional fiber of 10 km length. 
This does not only increase the demand for space, but 
also significantly decreases the maximal pulse repetition 
rate as well as the transmission. In the next section we 
discuss the realization of an experiment using detectors 
with significantly lower afterpulsing permitting to use 
smaller deadtimes and hence making smaller loops and 
higher pulse repetition rates applicable . 



B. Rapid gating detectors 

Rapid gating detectors [13, [IB [13 very short gates 
in order to limit the total avalanche current. This drasti- 
cally decreases afterpulsing and permits high gate rates. 
Practical deadtimes of about 10 ns allow to realize a time- 
multiplexer with fiber loops in the meter range. We im- 
plement 4 loops of 2m, 4m, 8m and 16 m generating 32 
bins in total. The overall transmission loss, averaged over 
all bins, is 1.35 dB, mainly due to the large amount of 
splices. 
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A sine wave generator (/ = 818 MHz, Agilent) is used 
as the central synchronization unit (FigH]). It applies 
gates to two Peltier cooled (-40°C) JDSU InGaAs APDs 
(effective gate width of 200 ps) and triggers the diode 
laser (1550 nm, 30 ps pulses, PicoQuant) with a heavily 
divided frequency. Details about the gating technique 
can be found in [3, [lH . In addition to these signals we 
generate an auxiliary signal at //8 = 102.25 MHz whose 
period of 9.78 ns corresponds exactly to the distance of 
adjacent bins. It is synchronized with the arrival of the 
bins and used to perform coincidence counting. By this 
means, detections that occur between bin arrivals are 
filtered out. The total length of the bin train is 150 
ns, which permits pulse repetition rates of up to 6 MHz. 
The number of coincidence detections for each laser pulse 
are registered on an oscilloscope (Lecroy 8600A). These 
numbers are recorded and analyzed later on. 

The dark count probability per coincidence gate, at a 
detection efficiency of 77 =16.5%, is equal to 1-10~^ for the 
first and 5 ■ 10^^ for the second detector. This efficiency 
value is close to the practical limit of our detectors and 
even a very small increase of the bias voltage leads to a 
significant growth of dark counts. With this setting we 
thus obtain a probability of 1 • 10~^ for a dark count to 
appear in a single-shot measurement (32 bins), which is 
negligible. 

The full characterization of afterpulsing is complicated 
[Tsj . Here we perform a relatively simple measure- 
ment to roughly estimate the impact of afterpulsing. We 
use a single detector and send pulses with /i =100 at a 
frequency of 512 kHz, which is well below the maximum 
of 6 MHz. Now we regard the number of detections in 
the 16 coincidence gates which are synchronized with the 
arriving bins and the number of detections in the imme- 
diately following 16 coincidence gates where no light is 
awaited. After 3000 pulses the sum of detections that 
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FIG. 3: Probability distributions p(nj^) of the number of 
detections n, when coherent pulses of average photon number 
fj, = 7, 14, 55 and 105 are sent onto our first time-multiplexing 
detector. 



FIG. 4: Timing schematic of rapid gating time-multiplexing 
detector. Central timing unit is a sine wave generator (/ = 
818 MHz) that is responsible for detector gating as well as for 
laser triggering (with a divided frequency (-^)). An auxiliary 
signal at f/8 is used to perform coincidence counting, working 
as an additional filter against noise. The coincidence signals 
are then registered in an oscilloscope, which is triggered at 
the laser repetition frequency (connection not drawn). 
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occurred in the first 16 gates is 12806, while it is equal to 
144 in the second 16 gates. This would thus add 1% of 
detections to the succeeding 16 coincidence gates in the 
limiting case of 6 MHz repetition rate, when signal pho- 
tons are possibly present in each coincidence gate. For 
most applications this contribution is insignificant. 

As done in the previous implementation, we character- 
ize the rapid gating time-multiplexing detector by send- 
ing weak coherent pulses with different mean number of 
photons /i. The acquisition using the oscilloscope takes 
a few minutes in order to obtain sufficient statistics. Re- 
sults for fj. =10, 50, 100 and 400 are shown in Fig\S\ 
The Monte Carlo simulation now includes an additional 
adaptation of the detection efficiency at high /it-values. 
An undershoot effect after a strong avalanche, caused by 
more than one photon, can lead to a non-discrimination 
of an avalanche in the succeeding bin. To obtain the op- 
timal fitting result the single photon detection efficiency 
has to be adapted from 16.5 % for low, to 14.5 % for high 
/i-values. Afterpulsing is not accounted for and no offset 
like in the conventional gating case (Sec lII A] ) is observed. 
With this adaptation, measurement results and simula- 
tion are in quite good agreement, showing in particular, 
that afterpulsing does not play an important role. 



III. ENERGY RESOLUTION 
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In the following discussion we will concentrate on the 
rapid gating implementation due to its higher efficiency 
and dynamic range. The same analysis can of course be 
performed in the same manner for the other implemen- 
tation. 

The complete set of conditional probabilities p{n\fi) 
would perfectly characterize our detector. However, 
even the measurement of closely spaced discretized /i G 
{1,2,3, Umax} would be very time consuming and 
not practical. Since our simulation agrees very well 
with the experimental data, we fill the detector re- 
sponse matrix M, where an element M^_„ := p{n\fj,), /i G 
{1, 2, 3, /imaa;}, with Calculated values between the 
measured supporting points =10, 50, 100, 200 and 400. 
The upper plot of FiglS] shows a graphical representation 
of M up to ^ = 150. 



A. Single-shot resolution 

Graphically it is quite evident that one can read M 
as well the other way round, meaning that we assume 
a certain number of detections n and want to infer an 
estimation of the average number of photons fj, of the 
initial pulse, i.e. the conditional probabilities p{fi\n) for 
a given n. 

Formally we infer the p{fj,\n) by using the Bayesian 
theorem and assuming a uniform distribution for p{p) in 
the range of characterization (/i = 0, 1, ■■■,Hmax)'- 
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FIG. 5: Probability distributions p(n|^) of the number of 
detections n, when coherent pulses of ^ = 10, 50, 100 and 400 
are sent on the rapid gating multiplexing detector. 
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Consider for example the distribution (bottom 
plot of Figini): the total probability for a > 100 to 
produce 1 detection is so small that plotting p{fJ.\l) for 
/^mos = 100 and fimax = 400 leads to the same result, 
thus is stable for a ^max > 100. In our case {^max = 400) 
the same argument holds for up to 15 detections. Single- 
shot events producing more than 15 detections need to 
be discarded, since the contributions from /i- values larger 
than 400 is not negligible. If one wants to use this detec- 
tor for pulses with higher average photon number, ^imax 
has to be adapted accordingly. We note that if respecting 
this rejecting rule, the assumption of a uniform distribu- 
tion of /i- values is the most conservative {— no informa- 
tion about incoming pulse, except /i < fimax)- 

FigEl (bottom) reveals that even with a single pulse 
measurement (number of detections between 1 and 15), 
the range of possible average photon numbers /i is quite 
restricted. For example in case of only one detection 
(p(/i-|l)), one can state on a 90 % confidence level, that 
the initial pulse had a ^ of 8 ±5^, corresponding to an en- 
ergy resolution of 4.2 attojoule (equivalent to 33 photons 
at 1550 nm). 

High resolution classical Joulemeters exhibit energy 
resolutions of down to 10 femtojoule, but are limited 
by pulse repetition rates of a few kHz High speed 

Joulemeters can process rates up to 100 kHz, but have 
energy resolutions on the nanojoule level. Our detector 
combines at the same time, high repetition rate of up to 6 
MHz and high single-shot energy resolution on the atto- 
joule level. We note that due to the short loop lengths of 
the rapid gating multiplexer (< 30 m), the mode disper- 
sion in a 50 fim core (NA=0.2), graded index multimode 
fiber amounts to less than 10 ps and can therefore be 
used to increase the amount of captured light. 



B. Multi-shot resolution 



FIG. 6: Top : Graphical representation of detector response 
matrix AI for coherent pulses from fi—Q to 150. Bottom : 
Distribution of /i-values for single-shot detection events from 
1 to 15 detections. 
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In the last step the assumption of uniformly distributed /i 
is used. The lower plot of FiglS]shows these distributions 
for different n from 1 to 15. 

The choice of iJLmax il^'max = 400 in our case) is of 
course arbritrary, but it is important to ensure that the 
distributions pifi\n) in EqU] are independent of fJ-max- 



Under the same assumptions as used before, we can 
obtain multi-shot distributions p(/i|ni,n2, ...,nfe), where 
ni, 77.2, ■■■,nk is a series of number of detections generated 
by coherent pulses with the same ^- value. Using the 
single-shot distributions p(/i|n) we calculate : 

k 

pifi\ni,...nk) ^ c-Y[pif^\ni) (2) 

2=1 

where c is a normalization factor ensuring that 
EJ^^g"" ■•■"-*:) — 1- The additional information 

from each shot drastically narrows down the /x-values 
that come into consideration. In FiglT] (top) we plot the 
evolution of the most probable /i- values when a series of 
5, 3, 4 and 5 detections is obtained. In the bottom plot 
the evolution of the relative error A^go%//i (90% con- 
fidence level), as a function of the number of shots for 
different coherent pulses (/i — 10, 50, 100), is shown. 

In addition to that, we plot the evolution of the rela- 
tive error achieved by a gated single-pixel detector (with- 
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FIG. 7: Top: Evolution of the most probable /i-values after a 
series of detection, here 5, 3, 4, 5. Bottom: Evolution of the 
relative error (assuming a 90 % width) with number of shots, 
when subjected to repetitive coherent pulses of /i = 10, 50 
and 100. As a comparison, we plot the relative error obtained 
by a single-pixel detector with variable attenuator to adjust 
a detection probability per gate of 50 % (see text for details). 



out multiplexer) with the same detection efficiency [t] = 
16.5%). The relation between the /i- value of the incident 
pulse and the detection probability per gate Pdet,gate is 
given by 



-1 
1] ■ a 



ln(l - Pdet.,gate) 



(3) 



The parameter a accounts for the case when an ad- 



ditional attenuator is used to avoid detector satura- 
tion. Pdet,gate Can be dynamically estimated via the ratio 
Ndet/Ngate of number of detections Ndet and the num- 
ber of appUed gates Ngate- Our results confirm the in- 
tuition that the relative error of /i decreases fastest for 
Pdet,gate about 50%. This could be achieved experimen- 
tally for pulses with /i > 4 using a variable attenuator. 
We neglect cases with /i < 4 in our discussion. The error 
on ^, i.e. A/igo%, is propagated from a 90% Poissonian 
error on Ndet- Due to the constant Pdet.gate the rela- 
tive error is independent of the strength of the incident 
pulses. The contrary happens in the time-multiplexing 
setup. The higher the /j,— value, the faster the relative 
error decreases. For example, at = 100 a relative er- 
ror of 0.1 is obtained after approximately 150 shots. The 
single-pixel detector needs roughly 4500 shots (not plot- 
ted) to achieve the same result. This implies a reduction 
of the number of shots by a factor of 30. 



IV. CONCLUSION 

We find that Peltier cooled InGaAs/InP APD, sen- 
sitive in the near-infrared wavelength regime (1100- 
1650 nm), are suitable for the implementation of time- 
multiplexing detectors. We realize two different versions. 
The first one is an easy-to-operate and easy-to-control 
version based on commercially available detectors. It has 
a large flexibility in choosing gate widths, which in return 
permits to process a large range of laser pulse widths (up 
to 100 ns). Due to afterpulsing one has to respect rel- 
atively long delays between adjacent bins which neces- 
sitates the use of rather long fiber loops (order of kilo- 
meters). The maximal pulse repetition rate is 22 kHz. 
The second implementation is based on rapid gating de- 
tectors, exhibiting very low afterpulsing. Deadtimes of 
10 ns allow the use of short fibers (a few meters) to re- 
alize the bin delays. 4 loops are used to create 32 bins, 
yielding a large response dynamic. The maximally at- 
tainable pulse repetition rate is 6 MHz and the effective 
gate width is fixed at 200 ps. Analysis of the detector re- 
sponse to coherent pulses makes it possible to calculate 
single-shot and multi-shot energy resolutions of pulses of 
unknown energy. We find that in single-shot events the 
resolution can be as small as 4.2 attojoule. 
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